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Foreword 



This document records three working studies from an ongoing 
investigation of models and methods for the prediction of the cumulative 
effect of weapons salvos. These papers are reproduced here in their entirety 
and in the chronological order in which they were produced. 

The first paper in the sequence, A comparison of an empirical rule for 
aggregating damage from a weapons salvo to a plausible model for the same 
purpose, resulted from an examination of an extant formula for estimating 
the expected proportion of damage to an area target from a weapons salvo. Its 
conclusion is that the formula gives optimistic, and in some instances 
impossible, results when compared to a “plausible" model for the effect of the 
salvo. 



The second paper. Damage aggregation for a weapons salvo by an 
empirical rule related to the Poisson approximation to the binomial, describes 
an alternative formula which is conservative when compared to the same 
“plausible" model. 

The third paper, A stochastic model for hit overlap in a weapons salvo 
directed against an area target that leads to a proportional mechanism for 
damage aggregation, describes a basic case of an emerging family of target 
configuration and weapons impact scenarios which lead to the damage 
aggregation mechanism of the “plausible" model. The “plausible" damage 
aggregation mechanism is renamed the proportional effects mechanism in 
this paper. 

The general setting for these studies is described in the first two 
paragraphs of the introduction to the first paper. The proportional effects 
(plausible) damage aggregation mechanism is derived in Section 3 of that 
paper. This information is recapitulated in the introduction to the 
subsequent papers. 
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A comparison of an empirical rule for aggregating damage 
from a weapons salvo to a plausible model for the same 
purpose 

1. Introduction 

The scenario considered here is that a salvo of n weapons is launched 
against a target. The number of weapons that hit the target is a random 
variable N with possible values 0, 1, ... , n. Possible damage to the target is 
measured as a percentage (or proportion) of the whole ranging from 0% to 
100%. The damage to a pristine target resulting from a single hit is a 
deterministic proportion d of the whole. The aggregate proportion of 
damage to the target from the salvo is a random variable D, the randomness 
in D resulting from the randomness in the number of hits N. 

The end objective is to predict the expected proportion of damage E(D) 
to the target resulting from the salvo. This prediction can involve modeling 
of the probability distribution of the number of hits N, and the way that 
damage aggregates as additional hits beyond the first are scored. For the 
immediate study, it is assumed that each weapon in the salvo hits 
independently with the same probability p, so that the number of hits N has 
a binomial probability distribution. This leaves the problem of how to 
aggregate the deterministic proportion of damage resulting from multiple 
hits on the target. 

The empirical rule for predicting E(D) is 
E(D) = 1 - (1 -d) E(N) 

where E(N) is the expected number of hits on the target. This rule ducks the 
issue of aggregating deterministic damage from multiple hits and to some 
extent ducks the issue of modeling the probability distribution for the random 
number of hits on the target. We will first examine the workings of the 
empirical rule in the case of a salvo of size one, and then proceed to a 
comparison of the empirical rule to the output of a plausible model for 
deterministic damage aggregation. 



2. For a salvo of size one 

For a salvo of size one, i.e. a single weapon, there is either a hit with 
probability p, or a miss with probability 1 -p. The consequences of the salvo 
are summarized in the following diagram. 
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Hit, N = 1, D = d 




Miss, N = 0, D = 0 



The expected number of hits is 

E(N) = 1 • p + 0 • (1 -p) = p 

and the expected proportion of damage is 

E(D) = d • p + 0 • (1-p) = pd 

The issue then becomes how in this case does the empirical rule 
E(D) = 1 - (1-df for predicting E(D) compare with the actual value 
E(D) = pd. The following arguments show that E(D) overestimates E(D ) 
for all values of p and d. 

We will hold p fixed and study the behavior of E(D) and E(D) as d 
varies from 0 to 1. First note that at d = 0 both E(D) and E(D) are 0. 

Then note that both E(D) and E(D) increase as d increases. This is evident 
by inspection of the two expressions, or from considering that the derivative 
of E(D) is 



_d_ 

dd 



E(D) 



_d_ 

dd 



{1 - (1 -d)P} 



and the derivative of E(D) is 



V 

(1 ~d)' V 



jjEiD) - -ijtfd) = p 

and that both derivatives are nonnegative. Also the derivative of E(D) is 
larger than the derivative of E(D), so that E(D) increases faster than £(D). 
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Since both expressions start at 0 for d = 0 and E(D) increases faster, it 
follows that E(D) is greater than E(D) for all values of d. 

It is worth noting that at d = 0, E(D) and E(D) have the same slope, 
i.e. 




d =0 




d =0 



P 



and that while the slope of E(D) continues at a constant p, the slope of E(D) 
grows to infinity as d approaches 1. The common slope of E(D) and E(D) 
at d = 0 suggests that E(D) can be a good approximation to E(D) for small 
values of d. This issue is studied numerically in the following four plots. 







3 



The empirical rule E(D ) becomes absurd as d approaches 1. At 
d = 1 it says that expected damage to the target is total regardless of the prob- 
ability that the single weapon scores a hit. The reality is that the expected 
proportion of damage can never exceed the hit probability. 

Wayne Hughes [1] makes the observation — "A single missile attack is 
most attractive tactically when the damage E(D) is thought to be high, that is 
when p and d together are high. We do not need to make a judgement as 
to what is "high" to conclude that the region of primary interest will be 
toward the right end of the plots, which is the region of greatest divergence 

between the empirical E(D) and the more plausible E(D) curves." 



3. A plausible model for damage aggregation 

The premise of the "plausible model” for damage aggregation is that if 
the proportion of a pristine target that is damaged by a single hit is d, then 
each additional hit damages the same proportion d of that part of the target 
not previously damaged. Thus if D(k) is the aggregate proportion of damage 
to a pristine target from exactly k hits, then 



D( 0) = 0 = 1- (1 -d)° 

D( 1) = d = l-(l-rf) 1 

D( 2) = D(l) + rf{l-D(l)} = d + dil-d) = 1 - (1 -d) 1 



D(n) = D(n-l) + d{l - D(n-l)} = 1 - (l-d) nA + d{\-d) nA 



= 1 - (l-rfXl-rf)"- 1 = 1 - (1 -d) n 



where n is the number of weapons in the salvo. 



It is worth noting that the incremental proportion of damage from the 
k^ hit is, for k = 1, ... , n, 

D(k)-D(k- 1) = {1 - (1 -d) k } - {l-(l-rf) M } = d{\-d) kA 
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4. For a salvo of size n 

For a salvo of size n in which each weapon hits independently with 
the same probability p, the random number of hits N has a binomial 
probability distribution, i.e. 

P[N=k] = ^ > jp k (l-p) n ~ k , k = 0, ... , n 
and E(N ) = up. 

Then the expected proportion of damage is 

n n 

E(D) = £ D(k) P[N=k ] = ^ {1 - il - d)k ) p [^=E] 

k= 0 k =0 



= 1 - ^ (l-d) k P[N=k] = 1 - X (1 ‘^ )A 

= 0 Jt=0 



= i - X (JnWo#* = 1 - S 

it =0 *=0 



= l - {(l-rf)p + (l-p))” = l-d-prf)” 

and the empirical rule for predicting E(D) is 

E(D) = l-(l-rf)"P 

As in the case of a salvo of size one, E(D) overestimates E(D ) for all 
values of p and d, since 

E(D) > E(D) » 1 - (l-rf)"P > 1 - (1 -pd) n <=> (1-prf)" > (l-rf)"P 
» (1-prf) > (1 -d)P <=> l-(l-d)P > pd 
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The inequality between E(D) and E(D) for a salvo of size n reduces to the 
inequality between E(D) and £(D) for a salvo of size one. 

As before, we will hold p fixed and study the behavior of E(D) and 
E(D) as d varies from 0 to 1. First note that at d = 0 both E(D) and E(D) 
are 0. Then note that both E(D) and E(D) increase as d increases. This is 
evident by inspection of the two expressions, or from considering that the 
derivative of E(D) is 

-fj E(D) = = np(\-d) n r> A 

and the derivative of E(D) is 

JjE(D) = jj{l-(l-pd) n ) = np(l-pd I”' 1 

and that both derivatives are nonnegative. At d = 0 both derivatives reduce 
to up, so that E(D) and £(D) start with the same slope. 

At d = 1, E(D) becomes 1 while E(D) is 1 - (l-p) n , the probability of 
scoring at least one hit. The probability of scoring at least one hit is a clear 
upper bound on E(D), even when total destruction from a hit is assured. 
Also, 



while 



dd 



E(D) 



d =1 



0 if n p > 1 
< 1 if n p = 1 
°° if tip < 1 




d =1 



np(l-p) n ~ l 



As the values of the derivative of E(D) at d = 1 suggest, there are 
three basic shapes that an E(D) curve can assume. These shapes are 
illustrated in the following plot. 
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That E(D) is concave in d for up > 1 and convex in d for np < 1 can be 
confirmed analytically by examining the second derivative of E(D) with 
respect to d. 




- np{np- l)(l-d)”P' 2 



< 0 if n p > 1 
>0 if np< 1 



On the other hand, E(D) is seen always to be concave in d by 
examination of its second derivative with respect to d. 



d 2 

— y E(D) = -M/>(w-l)(l-prf)”- 2 <0 
d d 

It is reasonable to suspect that E(D) will approximate E(D) poorly for large 
values of d when E(D) is convex and E(D) is concave, i.e when the curves 
break oppositely from their initial common slope. 
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The relative behavior of E(D) and E(D) in selected cases is shown in 
the following plots. 
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Damage aggregation lor a weapons salvo by an empirical 
rule related to the Poisson approximation to the binomial 

1. Introduction 

This working paper follows a previous working paper "A comparison of 
an empirical rule for aggregating damage from a weapons salvo to a plausible 
model for the same purpose" [1]. The material in the next three paragraphs of 
this introduction is summarized from the previous paper and is discussed in 
greater detail there. 

The scenario considered is that a salvo of n weapons is launched against a 
target. The number of weapons that hit the target is a random variable N 
with possible values 0, 1, ... , n. Possible damage to the target is measured as a 
percentage (or proportion) of the whole ranging from 0% to 100%. The 
damage to a pristine target resulting from a single hit is a deterministic 
proportion d of the whole. The aggregate proportion of damage to the target 
from the salvo is a random variable D, the randomness in D resulting from 
the randomness in the number of hits N. 

The premise of the model for damage aggregation is that if the proportion 
of a pristine target that is damaged by a single hit is d, then each additional hit 
damages the same proportion d of that part of the target not previously 
damaged. Thus if D(k) is the aggregate proportion of damage to a pristine 
target from exactly k hits, then 

D(k) = l - (l-d) k , k = 0 , ... ,n 

The objective is to predict the expected proportion of damage E(D) to the 
target resulting from the salvo. This prediction can involve modeling of the 
probability distribution of the number of hits N. If it is assumed that each 
weapon in the salvo hits independently with the same probability p, so that 
the number of hits N has a binomial probability distribution, then 

E(D) = 1 - (1 -pd) n 

The purpose of this paper is to show that when N has a binomial 
distribution so that the expected number of hits is E(N) = np, the empirical 
rule 



E(D) = i - e -EW)d 
is a conservative approximation to E(D). 
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Like the empirical rule considered in the previous paper, which is an 
optimistic approximation to E(D), this rule depends only on E(N). Since the 
exact formula for E(D) is simple in the binomial case, much of the interest in 
the empirical rules will focus on their behavior when the distribution of hits 
is not binomial. 



2. Derivation of E{D) 

The familiar Poisson approximation to the binomial distribution arises in 
the case that the number of trials (salvo size) n approaches infinity, while 
the probability of success (hitting the target) p on a single trial approaches 
zero in such a way that the expected number of hits maintains the constant 
value 



Then the binomial probability of exactly k hits approaches the Poisson 
probability 



If the number of hits N is assumed to have a Poisson distribution with 
parameter E(N) = A, then the expected damage to the target turns out to be 



A 



E(N) = np = A 



yk 

P[N=k] = — e‘ A 
k\ 



k = 0,1,... 



oo 



oo 






E(D) 



A 





* = 0 



*=0 




■j _ Q-dk — | _ Q-E(N)d 
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3. Properties of E(D) 

A 

That E(D) is a conservative approximation to E(D) when the true 
distribution of N is binomial with E(N) = np = X follows from the fact that 

(1 ~^d) n converges to e~ d ^ from below, so that E(D) = \-{\- ^d) n 
converges to E (D) = 1 - e~ d ^ from above, or from the sequence of 
analytical comparisons 

E(D) > E (D) « 1 - (1 -pd) n > 1 - e~ n P d 

<=> e' n P d > (l -pd) n <=> e'P d > 1-pd 

The final comparison is an application of a standard inequality about the first 
two terms of the power series expansion of e' x . 

As in the previous paper, we can hold n and p fixed and study the 
behavior of E(D) as d varies from 0 to 1. At d = 0, E (D) = 0. The 

A 

derivative of E (D) is 

— E(D) = np e~ n P d 
dd y 

A 

This derivative reduces to up at d = 0, so that E(D), E(D) / and the empirical 
rule considered in the previous paper, all begin at zero with the same slope. 

A 

It is clear by inspection that E (D) is increasing and concave. 

A 

The relative behavior of E(D) and E(D) in selected cases is shown in the 
following plots. 
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z.A stochastic model for hit overlap in a weapons salvo 
directed against an area target that leads to a proportional 
mechanism for damage aggregation 



1. Introduction 

This working paper follows two previous working papers "A comparison 
of an empirical rule for aggregating damage from a weapons salvo to a 
plausible model for the same purpose" [1] and "Damage aggregation for a 
weapons salvo by an empirical rule related to the Poisson approximation to 
the binomial" [2]. The material in the next three paragraphs of this 
introduction is in part summarized from [1] and is discussed in greater detail 
there. 

The scenario that has been considered is that a salvo of n weapons is 
launched against a target. The number of weapons that hit the target is a 
random variable N with possible values 0, 1, ... , n. Possible damage to the 
target is measured as a percentage (or proportion) of the whole ranging from 
0% to 100%. The damage to a pristine target resulting from a single hit is a 
deterministic proportion d of the whole. The aggregate proportion of 
damage to the target from the salvo is a random variable D, the randomness 
in D resulting from the randomness in the number of hits N. 

The premise of the model for damage aggregation has been that if the 
proportion of a pristine target that is damaged by a single hit is d, then each 
additional hit damages the same proportion d of that part of the target not 
previously damaged. Thus if D(k) is the aggregate proportion of damage to a 
pristine target from exactly k hits, then 



This mechanism for aggregating the cumulative effect of hits has been 
referred to as plausible in the previous papers. A more descriptive 
terminology would be to call it a proportional effects mechanism. 

The objective has been to predict the expected proportion of damage E(D ) 
to the target resulting from the salvo. Since 



D(k) = 1 - (1 -d) k , k = 0,...,n 
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where P[N=k] is the probability of exactly k hits from the salvo, this 
prediction involves modeling the probability distribution of the number of 
hits N. So far it has been assumed that each weapon in the salvo hits 
independently with the same probability p, so that the number of hits N has 
a binomial probability distribution. The impact of this assumption on E(D) is 
discussed in [1]. 

Effects mechanisms can be derived from assumptions about the 
geometry of the target, the coverage of the weapon, and the probabilities of 
hitting locations within the target area. Then the proportion of damage from 
k hits becomes a random variable A (k), and with D(k) = E{A(fc)} 

n n 

E(D) = E{A(k)}P[N=k) = £ D(k)P[N=k] 

k =0 k = 0 

Viewing the model at this deeper level of detail can provide a better picture of 
its applicability. 

The purpose of this working paper is to present one scenario for 
weapons overlap on an area target which leads to a proportional effects 
mechanism in the sense that 

D(Jfc) = E{A(fr)} = 1 - (1 -d) k , k = 0,...,n 

for a pertinent value of d and to examine the probability distributions for 
A (k) that it implies. Probability distributions for A (k) can be combined with 
probability distributions for N to obtain deeper rooted probability 
distributions for D, the proportion of damage to the target. Given a 
probability distribution for D, one can set a damage threshold sufficient to 
meet a tactical goal and predict the probability that the salvo will damage the 
target in the sense that D achieves the established threshold. 



2. A cellular area target scenario 

Suppose that an area target is divided into m disjoint cells. Each cell 
represents a portion of the target which would be damaged by a single weapon 
which impacts on that cell, so that m nonoverlapping hits (one on each cell) 
would exactly suffice to totally damage the target. 
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However suppose that each hit from the salvo impacts a randomly 
chosen cell within the area independently of the cells impacted by the other 
hits. Then if k hits impact j different cells, the proportion of the target 
which is damaged is 



Att) = 1 

If there are no hits ( k = 0), then A(0) = 0 and D(0) = E{A(0)} = 0. If 
there is only one hit (k = 1), then A(l) = 1/m and D(l) = E{A(1)} = 1/m. Thus 

D(0) = 1 - (l-d)° = 0 
D(l) = 1 -(1-rf) 1 = d 

for d = 1/m. These formulas for D(0) and D(l) hold regardless of the 
number of cells in the target. 

3. The single cell target 

Suppose that a target has only one cell (w = 1). Then the target is 
totally damaged by the first hit, and each subsequent hit has no additional 
effect. Thus Aft) = 1/m = 1, and D(Jt) = E(A (k)) = 1 = 1 - (1 -d) k 
for d = 1/m = 1, k = 1, 2, ... . Thus a single cell model leads, in a trivial way, to 
a proportional effects mechanism. 



17 



4. The two cell target 



Suppose that a target has two cells (m = 2). Then the possible 
proportions of the target that can be damaged after a sequence of hits are 
Mm =1/2 and 2/m =1. The aggregation of damage from successive hits can be 
represented by the simple transition diagram below 




1st hit 2nd hit 



3rd hit 4th hit 



where the possible proportions of damage after each hit are indicated by the 

corresponding values of A (k), the arrows indicate the possible transitions 
from each damage state to subsequent damage states, and the probabilities of 
such transitions are shown by labels on the arrows. 

The probability that the target is in damage state A (k) = 1 or A(k) = 1/2 
after k hits can be computed by multiplying transition probabilities together 
along each path leading to the state and then adding up the products. It 
follows that 

P[A(1) = 1/2] = 1 



so that 



D(l) = E{A(1)} = (1/2) P[A(1) = 1/2] = 1/2 = 1 -{1-0/2)} 1 
as previously observed. Continuing 

P[ A(2) = 1/2] = 1/2 
P[A(2) = 1] = 1/2 

so that 



D( 2) = E{A(2)} = (1/2) P[A(2) = 1/2] + (1) P[A(2) = 1] 
= (1/2) (1/2) + (1) (1/2) = 3/4 = l-{l-(l/2)} 2 



18 



Similarly 



P[A(3) = 1/2] = 1/4 
P[A(3) = 1] = 3/4 



so that 



D(3) = E(A(3)} = (1/2) P[A(3) = 1/2] + (1)P[A(3) = 1] 
= ( 1 / 2 ) ( 1 / 4 ) + ( 1 ) ( 3 / 4 ) = 7/8 = l -{ l -( l / 2)} 3 



In general 

P[A(A) = 1 / 2 ] = (1/2)** 1 
P[A(k) = 1] = 1 -(1/2) m 

so that 

D(k) = E{A(k)) = (1/2) P[A(k) = 1/2] + (1) P[A(k) = 1] 
= ( 1 / 2 ) { l / 2) t ' 1 + ( 1 ) ( l - d / 2))^ 1 
= 1 -( 1 / 2 )^' = 1 -{ l -( l / 2 )}^' 

Thus 

D(k ) = 1-d-d)* k = 0,1,... 



for d = l / m = 1/2 showing that the two cell model leads to a proportional 
effects mechanism. 

In the two cell case the probability distribution for Aik) has been 
particularly easy to compute. That distribution is tabulated below for the first 
few values of k. 



k 


0 


1 


2 


3 


4 


5 6 


7 


8 


P[A(k) = \/2) 


0 


1 


1/2 


1/4 


1/8 


1/16 1/32 


1/64 


1/128 


P[A(k) = 1] 


0 


0 


1/2 


3/4 


7/8 


15/16 31/32 


63/64 


127/128 
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5 . The m cell target 



Now consider a target with an arbitrary number m of cells, for the 
The transition diagram for the aggregation of damage from successive hits 
begins as shown below. The probability of an upward transition from any 
total damage state A (k) = m/m = 1 is zero, so that all upward transitions 
beyond total damage become impossible. 




1st hit 2nd hit 3rd hit 4th hit 



Again multiplying the transition probabilities along paths leading to a 
damage state and adding the products, it follows that 

P[A(l) = i] = 1 

so that 

D(l) = E(A<1) = JjP[A(l) = l] = X = l-(l-l) 1 
as previously observed. 
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Continuing 



P[A(2) = i] = i 
P[A(2) = i] = 



so that 

D(2) = £IA<2> = ±PlM2) = ±, + 2.p M 2) = il = i(l) + ^(a=i) 



Further 

P[A(3) = ±) = (i) 2 = -L 

P(A(3) = i] = (^i)(|) + (l)(^i) = 2<^i> 

««3)-il - (^)(^) = ( "'- 1 ,’ ( f- 2) 

so that 

D(3) = E{A( 3)} = ±P[A(3) = ±] + l f P[A(3) = £) + ^P[A(3) = ±] 



1 


r n 


,2. 


r 3 (m - 1)^ 


,3. 




_ 3m~-3m + l 


m 


l” j2 J 


+ m 


l m2 J 


+ m 


v nj2 ) 


m* 



The preceding computations of D(l), D(2), and D( 3) indicate that the 
damage aggregation mechanism for the m cell target is proportional until 
d = ]/m. It remains to confirm the proportionality of D(k) for an arbitrary k. 
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The general transition diagram, going from k hits to k + 1 hits, 
appears below 




A(Jt+l) = m/m 



MM) = (m-1 )/m 




A(Jt+l) = j/m 



MM) = (/-D/m 



A(fc) = 3/m 



A(fc) = 2/m 



A(fc) = 1 /m 




It follows that 

P[A(J; + l) = i] = ±p[A(k) = ±) 

and 

P[A(* + 1) = 1] = JjjP[A(k) = ±) + -- - Ji - - P[A(fc) = -£] ; = 2 m 
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Using the preceding k to k + 1 transition equations 



m m 

DU' + l) = E{A(fc + l)} = ^±p[A(k + l) = l] = ± P [ A (k + l) = ±] + ^±p[A(k + l) 

H j = 2 

m . - . . 

= W ■ 7H p [A«) = PtA(Jr) = -i] + P[A(« = £] 

,=2 

m . , m . 

;=1 /=2 



7=2 

m-1 



m . . m-1 . 

y=i ;=> 

m 



m m. 

/=1 ;=1 

m 



t 



m m . 

= Z ^ 4 p Wfc) = il + ^ 

7=1 7=1 



since the term indicated by an arrow is zero. Continuing 



m 



m 



Dflt + l) = 



XFsW)4)*X^ ! V p [«)=il 



H 

m 



M 



M 






m 



-i 



M 



i m-l j . 1 1 d [ a ( m - / 1 
| w m + m Lm/C ' mJ 



m 



rn 



m- 



H 



) 



(=1 



= -^ElAtt)! + Jf-1 = ^=J-D(« +■ 1 



m 



m 



m 



The recursive relation 



D(k + 1) = Z^lD(k) + ± 



which it turns out holds for k = 0,1, ..., permits the verification of 
proportionality, since if 



D(k) = 1- 




then 



D(k + 1) = 




and we have already verified proportionality for k = 1,2, and 3. This 
completes the demonstration that the damage aggregation mechanism for tl 
m cell target is proportional with d = 1/m. 

6. Damage distributions for the m cell target 

The k hits to k+\ hits transition equations provide a means for the 
calculation of probability distributions for A (k), k = 1,2,.... For this purpose 
these equations are perhaps more conveniently written as 



The numerical table which follows for A (k) when m = 3 was obtained by 
their recursive application. 



P[A(k + l) = i] = J f P[A(fc) = i] 




V J 



J 
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k 


i 


2 


3 


4 5 


6 


7 


8 


9 


10 


P[A(k) = 1/3] 


i 


1/3 


1/9 


1/27 1/81 


1/243 


1/729 


1/2187 


1/6561 


1/19683 


P[A{k) = 2/3] 


0 


2/3 


6/9 


14/27 30/81 


62/243 


126/729 


254/2187 


510/6561 


1022/19683 


II 

< 

a. 


0 


0 


2/9 


12/27 50/81 


180/243 180/729 


1932/2187 


6050/6561 


18660/19683 



Distributions for A(l), A(2), and A(3) were obtained in algebraic form 
in Section 5. They are repeated in the table which follows. The distribution 
for A(4) was obtained from the k = 3 to k = 4 transition equations. 




In the preceding tables it is understood that P[A(0) = 0] = 1 and 
P[A(k) = 0] = 0, fc = 1, 2, ... . 



Recursive application of the transition equations computes that 
P[A(k) = j/m]= 0 whenever j is greater than m. The formulas in the algebraic 
table reduce to 0 whenever j is greater than m. 

A possible algorithm for automating the computation of distributions 
for the A(k)'s is discussed in the appendix to this paper. 



Appendix - A posssible algorithm for computing damage probability 

distributions. 



The purpose of this appendix is to describe an approach to generating 
probability distributions for the A(J:)'s (the proportions of damage to an m 
celled target resulting from k hits) which has some potential for efficiency 
and the control of numerical error. 

The following expressions for PJ^-] = P[A(/c) = -£] were derived from 
the recursion equations at the beginning of Section 6, with the initial 

conditions PJ-^] = 1 and P,^] =0,j= 2, 3, ... , using the symbolic algebra 
software Theorist [3]. 






p 2 ta- 



rn - 1 



m 



P 3 fir] = 



__L 



P 3 [^=3 



m 



m - 1 



p 3 m = 



[m - 2][m - 1] 



m2 



P4(Jr]= J T P4ia = 7^— 1 
m 6 m 3 



[m - 2][w - 1] [m-3][m -2][m - 1] 

PaIM = 6 P 4 is - ] = : 



m 



m 



P 5 [Jt) = -L 



Ps(s) = 15—— P 5 a=25 lm - 211m - 11 



m 



m 



m 



p, [q _ jq [0! ~ [m ~ ~ ]} _ [rn - 4] [m -3][m -2][m-l] 



m 



m 



P 6 (^ = 65 



P6^]=- L -P6[H=31^— !- 
m 5 m 5 

[m - 3][m - 2][m - 1] 



Ps^ = 90 



[m - 2] [m - 1] 



m 



Pem = 15 



[m - 4] [m - 3] [m - 2] [m - 1] 



m 



m 



P 6 m = 



[in - 5][m - 4][w - 3][m - 2][m - 1] 



m 
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Each Pfc[£] has the structure 



PkU r] = c k [j] 



(m-l)(m-2) 



+ ^ 

— L.i 



rn 



k-1 



O'] 



ml 






where C k [j] is a coefficient depending only on k and j, a fact that can be 
confirmed by inspecting the m celled transition diagram at the beginning of 
Section 5. 

It appears that the generation of the P k [i]'s can be approached by 



seperately computing the terms — and the coefficients 

nr 1 

C k [j] using the modified version of the transition equations 

Q +1 [l] = C*[l] 

Q +1 [;1 = jC k [j] + C k [j-l] 



with the initial conditions CJl] = 1 and C } [j] = 0, j = 2 , 3 , ... . 



The following table of coefficients C k [j] was derived from the 
preceding coefficient transition equations. 



k 1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


c*[i] 


1 


1 


1 


1 


1 


1 


1 


1 


1 


Q[2] 


1 


3 


7 


15 


31 


63 


127 


255 


511 


Q-[3] 


0 


1 


6 


25 


90 


301 


966 


3025 


9330 


CJ4] 


0 


0 


1 


10 


65 


350 


1701 


7770 


34105 


Q[5] 


0 


0 


0 


1 


15 


140 


1050 


6951 


42525 


C,[6] 


0 


0 


0 


0 


1 


21 


266 


2646 


22827 


C k [7 ) 


0 


0 


0 


0 


0 


1 


28 


462 


5880 


Q[8] 


0 


0 


0 


0 


0 


0 


1 


36 


750 


C k l9] 


0 


0 


0 


0 


0 


0 


0 


1 


45 


CJ10] 


0 


0 


0 


0 


0 


0 


0 


0 


1 



27 



References 



[1] J. D. Esary. A comparison of an empirical rule for aggregating damage 
from a weapons salvo to a plausible model for the same purpose. Working 
Paper on Damage Aggregation, Naval Postgraduate School, April 1989. 

[2] J. D. Esary. Damage aggregation for a weapons salvo by an empirical rule 
related to the Poisson approximation to the binomial. Working Paper on 
Damage Aggregation, Naval Postgraduate School, April 1989. 

[3] Theorist. Presience Corporation, San Francisco, 1989. 



28 



INITIAL DISTRIBUTION LIST 



1. Library, (Code 0142) 2 

Naval Postgraduate School 

Monterey, CA 93943-5000 

2. Defense Technical Information Center 2 

Cameron Station 

Alexandria, VA 22314 

3. Office of Research Administration 1 

Code 01 2A 

Naval Postgraduate School 
Monterey, CA 93943-5000 

4. CDR W. A. Dewey 1 

Technical Operations Department (Code 80) 

Naval Undersea Warfare Engineering Center 
Keyport, WA 98345-0580 

5. James D. Esary 40 

Department of Operations Research 

Naval Postgraduate School 
Monterey, CA 93943-5000 

6. Wayne P. Hughes 1 

Department of Operations Research 

Naval Postgraduate School 
Monterey, CA 93943-5000 

7. Judith H. Lind 1 

Department of Operations Research 

Naval Postgraduate School 
Monterev, CA 93943-5000 

8. R. Neagle Forrest 1 

Department of Operations Research 

Naval Postgraduate School 
Monterey, CA 93943-5000 



9. Rodney L. Mash 1 

Proof, Test and Evaluation Department (Code 50) 

Naval Undersea Warfare Engineering Center 
Keyport, WA 98345-0580 

10. Gordon Nakagawa 1 

Department of Operations Research 

Naval Postgraduate School 
Monterey, CA 93943-5000 

11. Samuel H. Parry 1 

Department of Operations Research 

Naval Postgraduate School 
Monterey, CA 93943-5000 

12. Peter Purdue 1 

Department of Operations Research 

Naval Postgraduate School 
Monterey, CA 93943-5000 

13. Rex N. Randolph 1 

Weapons Effectiveness Branch (Code 3517) 

Naval Weapons Center 
China Lake, CA 93555 

14. LCDR Roger Stemp 1 

Operations Analysis Curricular Office (Code 30) 

Naval Postgraduate School 
Monterey, CA 93943-5000 

15. William G. Tonkin 1 

Target Vulnerability Office (Code 35104) 

c/o Code 3517 
Naval Weapons Center 
China Lake, CA 93555 

16. U. S. Army AMSAA 1 

Attn: AMXSY-GA (Will Brooks) 

Aberdeen Proving Grounds, MD 21005 



